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ABSTRACT 

In this paper we study the determinants of starless core temperatures in the Perseus molecular 
cloud. We use NIPj (1,1) and (2,2) observations to derive core temperatures (Tki n ) and data from 
the COMPLETE Survey of Star Forming Regions and the c2d Spitzer Legacy Survey for observations 
of the other core and molecular cloud properties. The kinetic temperature distribution probed by 
NH3 is in the fairly narrow range of ~ 9 - 15 K. We find that cores within the clusters IC348 and 
NGC1333 are significantly warmer than "field" starless cores, and Tkin is higher within regions of 
larger extinction-derived column density. Starless cores in the field are warmer when they are closer 
to class O/I protostars, but this effect is not seen for those cores in clusters. For field starless cores, 
Tkin is higher in regions in which the 13 CO linewidth and the 1.1mm flux from the core are larger, 
and Tkin is lower when the the peak column density within the core and average volume density of 
the core are larger. There is no correlation between Tkin and 13 CO linewidth, 1.1mm flux, density or 
peak column density for those cores in clusters. The temperature of the cloud material along the line 
of sight to the core, as measured by CO or far-infrared emission from dust, is positively correlated 
with core temperature when considering the collection of cores in the field and in clusters, but this 
effect is not apparent when the two subsamples of cores are considered separately. 
Subject headings: stars: formation — ISM: molecules 



1. INTRODUCTION 

Starless cores are the link between molecular clouds 
and Class protostars. They are smaller and denser 
(r ~ 0.1 pc, n > 10 4 cm~ 3 ) than, and kinematically 
distinct from, the molecular clouds in which they are 
found ([Goodman et alj I1998D . Because starless cores 
have no source of internal heating, the temperature of 
a starless core is determined by external illumination 
and self-shielding. The temperatures of starless cores 
are set by the interstellar radiation field, and in some 
cases luminous nearby stars, attenuated by the molec- 
ular cloud. The temperature profile within a starless 
core is determined by the heating of its outer layers and 
the internal extinction, making the centers of starless 
cores colder than their outer layers (Evans et al.l 120011 : 
IStamatellos et al.ll2007l IZucconi et al ll2001h . 

The temperatures of starless cores can be inferred from 
thermal emission from dust or from molecular lines. De- 
riving core temperatures from continuum emission is 
complicated by uncertainties in the dust emission spec- 
trum. For instance, in cold and dense regions dust grains 
are expected to coagulate and become covered in icy 
mantl es, creating variation in the emissivity spectral in- 
dex (|Ossenkopf fc Herm ing 1994). In addition, below a 
density of ~ 3 x 10 4 the d ust temperature is not cou- 
pled to the gas temperature (| Gol dsmith 2 00lHGalli et al.l 
2002), which is the more important quantity because the 
gas makes up ^99% of the mass. One of the best tracers 
of gas temperature is the rotation- inversion lines of NH3. 
Ammonia is particularly well suited for this type of study 
because 1) NH3 does not seem to deplete onto dust grains 
in starless cores, and in fact may have an enhanced abun- 
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dance where CO is depleted ([Tafalla et al.ll2002D 2) NH 3 
has hyperfine splitting of its inversion lines, making fits to 
the temperature, velocity and optical depth more certain 
3) the NH3 (1,1) and (2,2) transitions are separated by a 
small difference in frequency, so it is possible to simulta- 
neously observe both sets of lines. In this paper we use 
amm onia (1,1) and (2,2) observations ((Rosolowsk v et all 
12008ft to derive the temperature of starless cores. We ex- 
amine whether these data support the several models of 
the relationship between core propert ies and the molec- 
ular cloud where they are found (e.g.. lEvans et al.l 120011 : 
iBethell et al.l 12004 IStamatellos et aLll2004| ) 

Perseus is one of the Gould's belt molecular clouds 
included in the Spitzer legacy project "Fr om Molecu- 
lar C ores to Planet-Forming Disks" (c2d) ([Evans et al.l 
120031 ) and the COM PLETE survey of star-forming re- 
gions ([Ridge et al.l l2006f) . at a distance of 250 pc 
([Cernis fc Straizvsl 120031 : iBelikov et all I2002T ) . Because 
it is so well observed, there is a wealth of data on the 
cores and the cloud that enables us to study the influ- 
ences of various internal and environmental factors on 
the temperature of 50 starless cores in Perseus. 

2. DATA 

2.1. Far-Infrared 

As described in ISchnee et al.l (|2008l ), we used far- 
infrared emission maps at 70, 100 and 160 /im to derive 
the dust temperature and column density of the Perseus 
molecular cloud. The 70 an d 160 11m maps wer e obtained 
as part of the c2d survey ([Rebull et al.l l2007h . and the 
100 /im map was ta ken from the IRIS reprocessing of th e 
IRAS all-sky maps ([Miville-Deschenes fc Lagachell2005h . 
The resolution of the emission-derived column density 
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(Nh) and dust temperature (Td) maps with which we 
will compare the NH3 temperature (Tkin) is 40", as set 
by the 160 /im data. 

The dust temperature is calculated from the ratio of 
the observed flux at 70 and 160 /im, taking into account 
the color correction at both wavelengths and the emis- 
sion from transiently heated very small dust grains at 
70 /an. The column density is derived from the flux at 
160 /im and the dust temperature, allowing for a vari- 
able dust emissivity and normalized to the extinction 
measured from NIR reddening. Maps of the emission- 
derived dust temperature and column density, and the 
proce dure used to make t he maps, are presented in de- 
tail in [Schnieitin (fUl). 

2.2. Near- Infrared 

Data from the final release of the Two Micron All 
Sky Survey (2MASS) point source catalog was used to 
create a large-scale extinctio n map of Perseus as part 
of the COMPLETE survey (jRidge et alj [20061 ). The 
NICER algorithm estimates the column density towards 
background stars by comparing th eir observed colors 
with the (assumed) instrinsic colors (|Lombardi fc Alvesl 
2001). Although we expect the extinction map to be 
our least-biased est i mate of the cloud co lumn density 
([Schnee et all 120061 : iGoodman et~aTl 12008ft . its 5' reso- 
lution is roughly ten times coarser than our other maps, 
which smooths out density gradients and makes it some- 
what harder to compare the cloud and core properties. 

2.3. NH 3 (1,1) and (2,2) 

We used observations of 193 dense cores and core can- 
didates in the Perseus Molecular Cloud using the 100-m 
Robert F. Byrd Green Bank Teles cope (GBT) as de- 
scribed in iRosolowskv et al.l (|2008t ). For each target, 
single-pointing (31" FWHM) observations were made 
of the NH 3 (1,1) and (2,2) lines with 0.024 km s" 1 
spectral resolution. The targets for the ammonia sur- 
vey were drawn, in order of precedence, from 1) the 
locations of millimeter cores i n the Bolocam survey of 
the region ([Enoch et al.l l2006h 2) the locations of sub- 
millimeter cores in the SCUBA survey of the region 
(iKirk et alJ l 2006h 3) sources in the literature survey of 
Uijina et all (|1999h . and 4) cold, hig h-column-den s ity ob - 
jects in the dust map produced bv lSchnee et al.l (|2008f h 
and 5) sources identified by eye that appear in both the 
Bolocam and SCUBA maps but were not included in the 
published catalogs. The starless cores discussed in this 
paper were drawn from from only the first two categories. 

The kinetic temperatures of the starless cores (Section 
[3] describes how we identify which cores are starless) are 
determined by a simultaneous fit to the (1,1) and (2,2) 
emission profiles, accounting for collisional deexcitation 
of the two states. The same fit also yields the width of 
the NH3 line, which is decomposed into its thermal and 
non-thermal components based on the temperature de- 
termination. Although there arc likely gas temperature 
and density gradients within each st arless core on bot h 
theoretical and observational grouds ([Crapsi et al.ll2007h . 
the fitting of the ammonia-derived temperature is essen- 
tially an average over the core properties. However, a 
single-temperature model with uniform excitation condi- 
tions is a reasonably good fit to nearly all of the starless 
core in this sample. The NH3 observations can be used 



to accurately derive the kinetic temperatures of starless 
cores with 8 K < Tkin < 40 K where the lower and upper 
limits are set by sensitivity and transition t hermalization 
respectively (TWalmslev fc Un gcrcchts 1983). Further de- 
tails abou t the observations and d ata reduction are de- 
scribed in IRosolowskv et ail (|2008f ). 

2.4. 12 ' 13 CO (1-0) 

Observations of the 12 CO (1-0) and 13 CO (1-0) transi- 
tions at 115 and 110 GHz were carried at the 14-m Five 
College Radio Astronomy Observatory (FCRAO) tele- 
scope as part of the COMPLETE survey. The data have 
0.07 km s _1 spectral resolution and 40" spatial resolu- 
tion. The d ata reduction and mapping techniques are 
described in Rid ge et alj (|2006fh The CO maps cover 
approximately the same area in Perseus as the c2d sur- 
vey, giving us a good estimate of the kinematics of the 
ambient^ less dense, material in the molecular cloud. 



iPineda et all (|2008l) use this data set to determine the 
gas temperature and column density in regions in which 
the 13 CO (1-0) transition is optically thin (the 12 CO line 
is assumed to be optically thick everywhere) and the gas 
is not subthermally excited. 

2.5. 1.1mm Continuum Maps 

As part o f a su rvey that is complemetary to c2d, 
lEnoch et~aT1 {2006) mapped 7.5 deg 2 of the Perseus 
molecular cloud at 1.1mm with Bolocam at the Caltech 
Submillimeter Observatory. The resolution is 31", which 
matches very well with our NH3 observations, and all 122 
l.lmm-detected cores were observed with the GBT. The 
Bolocam observations are most sensitive to dense cores 
(n > 2 x 10 4 cm~ 3 ) with density contrasts at least 30 
times greater than the background cloud, which matches 
the description of a starless core embedded in a molecu- 
lar cloud quite well ([Bergin fc Ta falla 2007) . In addition 
to using the 1.1mm data for target selection, we use the 
1.1mm derived fluxes and sizes to compare the T^n with 
core properties such as mass, column density and volume 
density. 

3. IDENTIFICATION OF STARLESS CORES 

To ensure that we only study the starless core pop- 
ulation in Perseus, we compare the positions of our 
NH3 observations with the coordinates of known pro- 



tostar s, as found throug h the a nalysis of Enoc h et al 
(|2008[) . Uorgensen et al l (l2007f) and lHatchell et al 



(|2007l i. Uorgensen et alJ (|2007l ) identify 49 deeply em- 



bedded Young Stellar Objects (YSOs) in Perseus us- 
ing Spitzer/IRAC color cuts, Spitzer/MIPS and 850 
jim SCUBA correlat ions, and SCUBA morphologies. 
lHatchell et all (j20f)l uses NIR/MIR (2MASS/IRAC) 
detections and/or detections of protostellar outflows to 
separate starless cores f rom p rotostars, and find 56 pro- 
tostars. lEnoch et alJ (|2008h identifies 55 protostellar 
cores out of the total 122 millimeter cores in Perseus 
using the spatial coincidence between 1.1mm continuum 
cores and sources identified in the c2d Del ivery Docu- 
ment ([Evans et al.ll2007t lHarvev et al1[200l as "YSOc" 
(YSO candidate), "red" or that are 70 //m point sources 
and are not classified as "Gale" (Galaxy candidate). In 
this paper we only consider those 50 ammonia pointings 
that are called "starless" by the various definitions of 
all three papers, or that are more than 30" from any 
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identified protostar in any of the three papers. None of 
the cores we classify as "starless" are identified as a low - 
luminosity embedded protostar by lDunham~e t al. (2008). 
We also require that the starless cores have solid NH3- 
derived temperatures and linewidths, meaning that T^m 
is greater than 2.73 K, the NH3 (2,2) detection is > 5cr 
and the model and observed spectra fit well by eye. Of 
the 193 observations in the NH3 survey, 137 of the point- 
ings are not coincident with an identified protostar, and 
50 of these satisfy the signal to noise and goodness of fit 
requirements. 

4. ANALYSIS 

Here we compare the ammonia-derived temperature 
with several properties of the starless cores and their en- 
vironment. We considered only those 50 starless cores 
with well-defined T\i n and linewidth, so we are able to 
derive their non-therm al linewidths. All w ere detected 
in the 1.1mm survey of lEnoch et alJ (|2006l ). so for each 
we have the total flux, peak flux (from which we derive 
the peak column density), mass (from the total flux and 
the NH3 temperature) and mean particle density. The 
FIR-derived dust temperature and column density are 
considered to be believable as long as they fall within 
the range 12 < T d < 20 and < A v < 24.5, which is 
true except in those regions where the 160 /jm map is sat- 
urated or corrupted by a nearby bright source. There are 
46 starless cores with believable FIR-derived Td and Ay , 
and all 50 starless cores fall within the area of our NIR- 
derived extinction map. The gaussian fitting algorithm 
failed to converge when calculating the 13 CO linewidth of 
Too for 6 cores. The internal properties of our sample 
of starless cores are listed in Table [TJ and the proper- 
ties of Perseus cloud along the line of sight to the cores 
are given in Table [2] Although the distance from each 
starless core to the neares t protostar found in the sur- 
veys of lEnoch et all ()2008f ): iJorgensen et alJ (|2007t ) and 
lHatchell et al.l (|2007l) is known, in our discussion of pos- 
sible correlations we consider only those with separations 
of 0.25 pc or smaller. 

4.1. Tldn 

The kinetic temperature distribution probed by NH3 
is quite narrow, from ~ 9 - 15 K, as seen in Fig. [TJ This 
may be due to the relatively low angular resolution of the 
GBT pointings, whose 31" beam corresponds to probing 
the inner 0.019 pc (3.9 xlO 3 AU) of each core. Models of 
externally heated prestellar cores show that they are only 
significantly colder than 9 K within the inner ~ 10 3 AU 
or warmer than 14 K outside a radius of ~ 10 4 AU, with 
variation dependin g on the strength of the ISRF inci- 
dent upon the core (|Evans et al.ll200lt IStamatellos et al.l 
l2004ft . We are unlikely to be biased by chemical effects, 
given that NH3 is not s een to deplete in starless cores 
(e.g., iTafalla et all 12002ft The critical density of NH 3 , 
~ 10 4 cm~ 3 (|Swadej 1 1989ft . is low compared the typical 
gas densities observed in the cores in our sample (10 4 
cm -3 < n < 10 6 cm -3 ), so we expect the NH3 observa- 
tions to sample the full volume of our cores. 

The range of temperatures seen in the Perseus starless 
cores is typic al of those seen in other molecular clouds. 
For instance, Kir k et all (|2007l ) derived dust tempera- 
tures of pre-stellar cores in Taurus and Ophiuchus and 
found that most cores fall within the 9-13K range. The 



peak of the T^m distribution shown in Figure[T]is between 
10-11K, which agrees well with the median temperature 
(11 K) derived from NH3 obs ervations of the 36 starless 
as observed and identified by iJiiina "eTaLl (fl999ft across 
many nearby star forming regions. T he lower and upper 
quartiles of the temperatures of the IJijina et alJ (1999) 
starless core sample (9.7 K and 14.0 K) are also a close 
match to that shown in Figure [1] 

Numerical simulations of starless cores have been used 
to quantify how their temperature profile is affected by 
variations in many parameters, such as the strength of 
the ISRF, the density profile, and the amount of ex- 
tinction provided by the molecular cloud in which the 
cores are embedded (e.g. lEvans et alj|2001t iGalli et al. 
2002t IStamatellos~fc Whit worthll2003t IStamatellos et al. 
2004ft . lEvans et al.l (|2001ft calculate that for Bonnor- 
Ebert spheres with central densities varying by about 
a factor of 100 (from 10 5 to 10 7 cm -3 ), the central tem- 
perature varies by a factor of ~2 (6.5 - 10 K), but the 
temperature at a radius of 10 3 5 AU (the radius of the 
GBT beam at the distance of Perseus) only varies by ~1 
K (9 - 10 K). A doubling of the strength of the ISRF 
results in a 1 K increase of the core temperature file over 
the ra nge of radii probed by the NH3 observations in this 
study (|Evans et al.ll2001ft . while a factor of 100 increase 
in the ISRF would result i n an increase in t he core tem- 
perature by a factor of ~2 (|Galli et al.ll2002l ) . Variations 
in the du st opacity also affect the temperature of starless 
core, but iGoncalves et alJ (|2004l ) show that this effect is 
quite small (< 1 K) for the range of dust properties that 
one expects in these cores. Given that the temperature 
of a prestellar core (at a radius of 10 3 5 AU) is not es- 
pecially sensitive to its central density, the dust opacity, 
the strength of the ISRF or the degree of extinction for 
the range of values that we expect in Perseus, the narrow 
width of the temperature distribution shown in Figure [1] 
and its peak at 10 K are roughly what one would expect. 

4.2. Tkin vs. Location 

The presence of local star formation appears to affect 
Tkin- Fig. [1] shows the temperature distribution of all 
the starless cores in Perseus (filled grey histogram), as 
well as that for starless cores in the clusters NGC1333 
and IC348 (open solid histogram) and field starless cores 
(open dashed histogram). The boundaries for NG C1333 
and IC348 are taken from I J0rgensen et alJ ((2007) . It is 
clear that those cores in clusters are warmer than those 
outside clusters. The probability that the temperature 
distributions of starless cores inside and outside clus- 
ters have the same mean is ~ 10~ 5 , as calculated by a 
Wilcoxon Rank-Sum Test, Student's T-statistic or two- 
sided Kolmogorov-Smirnov test. 

There are reasons to expect that starless cores in clus- 
ters would be warmer than those outside clusters. For 
instance, a cluster of protostars may directly heat their 
environment. In addition, protostars inject turbulence 
into the surrounding medium, making it more porous 
and less effective a s a shield from the ISRF. In addition, 
IJijina et"al . (1999) find that dense cores that are not as- 
sociated with IRAS sources are warmer when found in- 
si de clusters than outside. 

IStamatellos et all (|2007l ) predict that the presence of 
a nearby Class protostar should not affect the temper- 
ature of a nearby starless core at all, and that a nearby 
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Class I protostar will only heat the outer layers of a core, 
and only by ~ 1 — 2K. We test this by comparing the 
NH 3 -derived temperatures of starless cores in Perseus to 
the distance t o the nearest Class and I protostars (as 
determined bv lEnoch et all (|2008l )). As shown in Fig. [2] 
and Table El we confirm that there is no significant cor- 
relation between Tkin and the distance to the nearest 
protostar for those starless cores in the clusters NGC1333 
and IC348. However, there is a significant negative corre- 
lation between Tkin and the distance to the nearest class 
and/or I protostar for those starless cores in the field 
(i.e. starless cores further from protostars are colder). It 
is not at all certain that this result contradicts the pre- 
diction of IStamatellos et all (|2007l ). because we cannot 
say that nearby protostars heat starless cores, just that 
warmer cores are found near protostars. 

4.3. Ikin vs. Cloud Temperature 

Here we compare the temperatures of starless cores 
with the temperature of the material in which it is em- 
bedded. We estimate the line-of-sight avera ged tem- 
perature of the dust (Td) (iSchnee et al.l l2008h and the 
gas (Too) ([Pineda et al.ll2008l ) in the Perseus molecular 
cloud and compare it with the temperature derived from 
the NH 3 pointings. The resolution of the Td and Tco 
maps are both ~40", which is roughly comparable to the 
31" resolution of the ammonia data. However, both Td 
and Tco are imperfect measures of cloud temperature. 
Because 12 CO is optically thick (and perhaps depleted) 
it is not sensitive to the colder and denser cloud material 
in which the starless cores are embedded. The dust tem- 
perature determination is emission- weighted, and at the 
FIR wavelengths (70 and 160 /im) of the Spitzer maps 
the emission is dominated by warm dust along the line 
of sight, and therefore is less sensitive to the colder and 
denser clumps in which cores form. 

Nevertheless, one might expect that the temperature of 
a starless core could be affected, or predicted, by the tem- 
perature of the medium in which is it embedded. Starless 
cores are, by definition, externally heated, so whatever is 
heating the cloud should also be heating the embedded 
cores. This analysis can be biased because the gas and 
dust derived temperatures are warmer within the clusters 
than outside, which is also true for Tkin, so any derived 
correlations must take this into account. Figure [2] and 
Tab. [3] show that although both Tco and Td are signif- 
icantly and positively correlated with Tkin when consid- 
ering the entire population of starless cores, this effect 
is entirely explained by the cluster/field bias, and within 
each sub-population there is no significant correlation. 

4.4. Tkin vs. Column Density 

The heating of starless cores is governed by the 
strength of the local interstellar radiation field (ISRF). 
Given that there are no O or B stars within the Perseus 
molecular cloud, the ISRF seen by the cores should de- 
pend largely on dust shielding but not on the distance to 
any special source of external radiation. Because we are 
unable to measure the volume-averaged shielding around 
each starless core, we use as a proxy the column density 
along the line of sight. Based on this simple analysis, one 
would predict that Tkin would be anti-correlated with 
column density. 



On the other hand, for starless cores that are in pres- 
sure equilibrium, the inward pressures due to the weight 
of the Perseus molecular cloud and the self-gravity of 
the cores will be counter balanced by the internal gas 
pressure from thermal and non-thermal motion. The ex- 
ternal pressure on the starless cores from the molecular 
cloud is proportion al to the square of the column den- 
sity ( Pcioud oc Ay) (jBertoldi fc McKedll992l lLada et all 
2008) and the internal gas pressure from thermal mo- 
tions is proportial to Tkin (^thermal oc (T 2 ^ oc Tkin)- From 
a balance of pressure analysis, one would expect a posi- 
tive correlation between Tkin and column density, though 
this would not necessarily be a strong effect because out- 
ward pressure is also supplied by non-thermal motions 
and magnetic fields. We calculate that for the starless 
cores in our sample the inward pressure of the weight of 
the molecular cloud is roughly a factor of 2 weaker than 
their self-gravity. 

Here we compare the temperatures of starless cores 
with the column density of molecular cloud material 
along the line of sight to each core. We estimate 
the cl oud column density using far-infrared dust emis- 
sion (ISchnee et al l l2008fi and near-infrared extinction 
(|Rideeet all 12006). Although extinction may be the 
most secu re estimator of colum n density in a molecu- 
lar cloud (Go odman et al.l 12008). the 40" resolution of 
the emission-derived column density map compares fa- 
vorably with the 5' resolution of the near-infrared ex- 
tinction map. As with the cloud temperature, there is 
an inherent bias in that the column density in the clus- 
ters is higher and the core temperatures are higher, so 
we consider the field and cluster starless cores separately. 
Figure [Hand Table [3] show that there is a significant pos- 
itive correlation between extinction-derived column den- 
sity and Tkin when considering just the starless cores in 
the field and just the starless cores in the clusters, as 
well as for the ficld+clustcr ensemble. There is no signif- 
icant trend between Tkin and the dust emission-derived 
column density for field or cluster cores. 

One should keep in mind that the emission-derived col- 
umn density derived from dust emission loses accuracy 
at values of Ay > a few, so it is not so surprising to get 
different res ults depending on h ow the column density 
is estimated (|Schnee et al.l l2006). The dust temperature 
derived from FIR maps is emission-weighted, and the 
column density derivation depends on the dust temper- 
ature, so the column density along lines of sight with 
significant variations in temperature and density is both 
systemically underestimated and subject to large scatter. 
This effect could wash out any true correlations between 
the column density of the cloud material and the tem- 
perature of the core embedded in it. 

In a numerical simulation of e mbedded starless cores, 
IStamatellos fc Whitworthl (|2003l ) show that the temper- 
ature profile of a core shielded by Ay = 5 material (a 
typical value for Ay seen towards the starless cores in 
Perseus) is not significantly different from a naked core, 
while deeply embedded cores (Ay = 20) are noticably 
colder than naked cores. In addition, we expect that the 
column density along the line of sight is a poor indicator 
of volumetric shielding due to the porosity of molecular 
clouds like Perseus. The visual impression given by two- 
dimensional maps and position-position- velocity cubes of 
nearby molecular clouds is one of non-relaxed structures 
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with many filaments, voids and shells, supporting this 
view. The observational evidence supports the idea that 
the shielding provided by higher Ay material in Perseus 
does not lead to lower temperatures of the embedded 
starless cores. The opposite effect is observed, which may 
be explained by the greater thermal support needed to 
balance the increased pressure on the cores by the weight 
of the molecular cloud. 



4.5. Tkin vs. ctnt and aco 

Clumpiness formed by turbulence in a molecular cloud 
is known to significantly affect the temperature of the 
cloud material, with higher turbulence leading to greater 
porosity and penetrat ion by the ISRF, a nd therefore 
higher temperatures dBethell et all l2004h . However, 
IStamatellos et all (|2007l ) find that the densest portions of 
the cloud (where starless cores are found) are primarily 
heated by long wavelength radiation that is not effec- 
tively shielded by the clumpy structure. It is therefore 
not obvious if one would expect starless core tempera- 
tures to be correlated with turbulence. 

Here wc examine whether turbulence affects the tem- 
perature of the starless cores in the Perseus molecular 
cloud. We use the non-thermal lincwidth of NH3 as a 
proxy for the turbulence within the cores, and the total 
linewidth of 13 CO as an estimate of the turbulence in the 
surrounding medium (the measured 13 CO linewidth is 
roughly a factor of f higher than the thermal linewidth). 
As before, we consider field and cluster starless cores sep- 
arately to avoid confusing a correlation between Tkin and 
turbulence with a correlation between Tkin and the field 
vs. cluster environment. 

Figure [3] and Tabic [3] show that there is no correla- 
tion between the non-thermal linewidth of NH3 and Tkin 
when considering the field and cluster popoluations of 
starless cores separately, but there is a significant posi- 
tive correlation for the ensemble. This implies that those 
starless cores in Perseus are more turbulent when found 
in clusters than in the field, but the turbulence itself docs 
not necessarily lead to higher temperatures. Our expec- 
tation of starless core geometry is that they are in some 
sense "relaxed" structures and therefore not porous, re- 
gardless of the level of internal turbulence. CO linewidth 
does not correlate with Tkin for those cores in clusters, 
but it is significantly positively correlated with Tkin for 
field starless cores, as shown in Figure [U One inter- 
pretation of of the field/cluster difference is that within 
clusters the CO linewidth is likely to be affected strongly 
by the ongoing star formation, while the CO linewidth in 
the field better traces the cloud geometry, so only in the 
field does a larger CO linewidth imply better penetration 
by the ISRF. 



4.6. Tki n vs. 1.1mm Continuum Derived Properties 

It might be supposed that the temperatures of starless 
cores are partly determined by the internal properties of 
the cores themselves, and not just by their environment. 
Here we compare Tkin with f.Imm Bolocam flux, the 
peak column density (expressed in Ay), the mass and 
the average volume density derived from the 1.1mm map 
and NH3-derived temperature of each core. 

The peak Ay is derived from the peak 1.1mm flux 



beam 



) of each core and the NH3-derived temperature 

gbeam 



iV(H 2 ) 



(1) 



Qbeamtl'H2 m HK v B v (Tki n ) 

with N(E 2 )/Ay = 0.94 x 10 21 cm 2 mag" 1 
(|Frerking et al.1 [T"982f ). Here flbeam is the beam solid 
angle, m# is the mass of hydrogen, «i.i mTO = 0.0114 
cm 2 g _1 is the dust opacit y per gram of gas interpolated 
from Table 1 column 5 of lOssenkop f fc HennineJ (|1994j ) 
for dust grains with thin ice mantles, B v is the Planck 
function, and /i# 2 = 2.8 is the mean molecular weight 
per H 2 molecule. A gas to dust mass ratio of 100 is 
included in K\,i mm . 

The core mass is derived from the total 1mm flux (S v ) 
and Tkin using the equation, 

d 2 <? 

M= j,t\ > ( 2 ) 

where k„ is the dust opacity, d is the distance to Perseus 
(250 pc). The density is derived from the core mass and 
the deconvolved average FWHM size (D) of the core (n = 
6M/nD 3 n p ), where [i v = 2.33 is the mean molecula r 
weight per particle, as explained in lEnoch et al.l (2006). 

When deriving the core mass, peak column density and 
average volume density from the 1.1mm continuum emis- 
sion, we use the gas temperature (Tkin) in the calcula- 
tion, under the assumption that the gas and dust temper- 
atures are coupled. Although this assumption is valid for 
the densities (n > 2 x 10 4 cm" 3 ) d erived for these cores 
(iColdsmithl [20011 : iCalli eTaD 120021) . we are introducing 
a potential bias in our analysis because Tkin and the 
derived mass are not independent. It would be prefer- 
able to use the Tdust rather than Tkin to determine the 
l.lmm-derived core properties, but this requires map- 
ping the dust emission from each core at several wave- 
lengths, sampling both the peak of the spectral energy 
distribution (SED) as well as the Rayleigh- Jeans tail to 
derive the dust temperature, column desnity and emis- 
sivity spectral index profiles. Existin g surveys of Perseus 
([Enoch et al.ll2006t iKirk et al J [20061) have only sampled 
the Rayleigh- Jeans portion of the SED, so we do not have 
accurate estimates of the dust temperatures within the 
starless cores, and we are forced to use the gas tempera- 
ture in our calculations. 

Note that some of the Bolocam cores, whose numbers 
are given in Tabic Q] (Pcr-Bolo 24, 75, 100 and 122), have 
been observed in multiple GBT pointings with each posi- 
tion corresponding to different sub-condensations within 
the core. Each of these pointings is assigned the same 
1.1mm total and peak flux and dimensions as that of the 
entire Bolocam core, and the peak Ay, mass and den- 
sity are all derived from the overall core Bolocam flux 
measurements and the Tkin derived for the sub-core. An 
alternative analysis would have been to break up each 
Bolocam core with more than one NH3 pointing into sub- 
cores and assign different fractions of the total 1.1mm 
flux to each portion, but the procedure to do the divi- 
sion would be somewhat arbitrary, just as the decision 
of which cores should be observed only once and which 
in multiple positions was not based on an absolute pre- 
scription. Since neither analysis is obviously the correct 
one, we decided to go with the simpler of the two and 
double-count the 1.1mm properties of the four Bolocam 
cores with more than one NH3 observation. 
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In Fig. [3] and Table [3] we compare Tkj n with the 1.1mm 
flux and properties derived from the lEnoch et al.l (2006) 
Bolocam map of Perseus. We find no correlation between 
Tkin and the 1.1mm flux, core mass, average volume den- 
sity or peak Ay for those cores in IC348 and NGC1333. 
On the other hand, for starless cores in the field there 
are significant negative correlations between Tkin and the 
peak Ay and average volume density, and a positive cor- 
relation between Tkin and the Bolocam flux. If we as- 
sume that starless cores are effectively self-shielded from 
the ISRF because they are not porous, then we would 
expect higher peak column densities and volume den- 
sities to correspond to higher volumetric shielding and 
lower internal temperatures. It is not obvious why the 
self-shielding argument does not hold for those starless 
cores in NGC1333 and IC348. The connection between 
1.1mm flux and Tkin may be simply that for a given mass 
a warmer core will produce more emission, but this ar- 
gument does not explain why the positive correlation is 
not seen in cluster starless cores. 

5. CONCLUSIONS 

In this paper we investigate the temperature distribu- 
tion of starless cores in Perseus. Although, the range 
of temperatures derived from NH3 (1,1) and (2,2) point- 
ing is quite narrow, from ~ 9 — 15 K, there are several 
interesting correlations and lack of correlations. 

The location of a starless core does affect its tempera- 
ture. We find that those cores inside the two main star- 
forming regions in Perseus (IC348 and NGC1333) are 
warmer than those not in clusters. Although there is no 
correlation between Tki n and the distance to the near- 
est class 0/1 protostar for those starless cores in clusters, 
starless cores in the field are colder the further away they 
are from protostars. 

Certain properties of the molecular cloud are related 
to the temperatures of the starless cores embedded in 
it. Tkin is higher when the extinction-derived column 
density of the molecular cloud is larger, despite the ad- 
ditional shielding from the ISRF that the higher column 
density implies, which may be due to the increased pres- 
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Some internal properties of starless cores are corre- 
lated with their temperature. For those starless cores 
in the field, total 1.1mm flux is higher in cores that are 
warmer, and peak column density and average volume 
density are higher in cores that are colder. A higher 
temperature naturally leads to more 1.1mm flux for a 
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TABLE 1 

Internal Properties of Starless Cores 



RA 1 


Dec 2 


Tkin 3 (tnt 4 1.1mm Flux 5 


Peak A v 6 


Mass 7 


Density 8 Number 9 


J2000 


J2000 


K km s — 1 Jy 


mag 


^^-solar 


cm -3 



3:25:07.8 


30:24:22 


9.2 





09 


0.1 


9.1 


0.3 


1.4E+05 


1 


3:25:09.7 


30:23:53 


9.2 





1 1 


0.1 


10.1 


0.1 


3.3E+05 


2 


3:25:10.1 


30:44:41 


10.0 





12 


0.4 


9.0 


0.9 


2.5E+04 


3 


3:25:17.1 


30:18:53 


10.3 





11 


1.2 


10.1 


2.8 


4.5E+04 


i 


3:25:26.9 


30:21:53 


9.1 





10 


1.2 


12.3 


3.4 


5.3E+04 


6 


3:25:46.1 


30:44:11 


11.2 





18 


0.9 


14.2 


1.8 


8.5E+04 


11 


3:25:48.8 


30:42:24 


9.1 





13 


0.5 


34.6 


1.3 


3.8E+05 


13 


3:25:50.6 


30:42:02 


9.0 





1 1 


0. 1 


29.6 


1.2 


3.2E+05 


1 1 


3:27:02.1 


30:15:08 


10.4 





13 


0.6 


7.8 


1.3 


4.2E+04 


19 


3:27:28.9 


30:15:04 


10.7 





10 


0.7 


7.3 


1.5 


5.1E+04 


20 


3:27:40.0 


30:12:13 


10.5 





17 


0.4 


18.6 


0.9 


1.6E+05 


23 


3:27:55.9 


30:06:18 


9.1 





05 


0.7 


18.9 


1.9 


1.9E+05 


24 


3:28:05.5 


30:06:19 


9.4 





11 


0.7 


17.7 


1.8 


1.8E+05 


24 


3:28:32.4 


31:04:43 


10.6 





13 


0.5 


7.6 


1.1 


7.1E+04 


26 


3:28:33.4 


30:19:35 


10.7 





12 


1.0 


10.0 


2.2 


3.7E+04 


27 


3:28:39.5 


31:18:35 


11.7 





17 


1.9 


21.2 


3.6 


2.3E+05 


31 


3:28:48.5 


31:16:03 


11.3 





10 


0.2 


7.6 


0.4 


1.8E+05 


35 


3:28:52.2 


31:18:08 


13.8 





17 


0.4 


6.8 


0.5 


1.4E+05 


37 


3:29:18.5 


31:25:13 


11.9 





12 


1.4 


18.1 


2.5 


1.1E+05 


53 


3:29:22.5 


31:36:24 


9.8 





12 


0.6 


7.6 


1.5 


4.6E+04 


56 


3:29:25.8 


31:28:17 


10.3 





11 


0.3 


18.8 


0.8 


6.3E+05 


58 


3:30:24.1 


30:27:39 


10.5 





07 


0. 1 


7.0 


0.8 


4.2E+04 


61 


3:30:45.6 


30:52:36 


10.5 





13 


0.7 


9.4 


1.5 


5.5E+04 


63 


3:30:50.5 


30:49:17 


9.8 





08 


0.2 


6.4 


0.5 


9.9E+04 


61 


3:32:26.9 


30:59:11 


10.5 





13 


0.9 


11.0 


1.9 


1.0E+05 


67 


3:32:44.1 


31:00:01 


10.2 





11 


2.0 


16.7 


4.6 


9.2E+04 


70 


3:32:51.3 


31:01:48 


11.2 





22 


1.1 


10.7 


2.3 


7.2E+04 


71 


3:32:57.0 


31:03:21 


10.1 





12 


0.9 


16.8 


2.2 


1.6E+05 


72 


3:33:02.0 


31:04:33 


10.0 





14 


0. 1 


18.8 


0.9 


1.8E+05 


74 


3:33:04.3 


31:04:57 


9.9 





18 


0.4 


19.2 


0.9 


1.8E+05 


75 


3:33:06.3 


31:06:26 


10.3 





15 


0.4 


17.7 


0.9 


1.6E+05 


75 


3:33:25.2 


31:05:35 


9.7 





10 


0.4 


11.1 


0.9 


1.5E+05 


82 


3:40:49.5 


31:48:35 


12.4 





08 


1.0 


7.6 


1.7 


3.5E+04 


89 


3:41:40.2 


31:58:05 


9.5 





07 


0.2 


7.5 


0.6 


8.7E+04 


92 


3:41:46.0 


31:57:22 


9.6 





07 


0.3 


10.2 


0.8 


1.0E+05 


91 


Q . A Q . A A fl 

o:4v5:44.U 


qo . no . c;o 
oz:U2:oz 


12.9 





13 


0.8 


11.7 


1.3 


1.5E+05 


100 


3:43:45.5 


32:01:44 


14.0 





10 


0.3 


5.3 


0,1 


8.0E+04 


101 


3:43:45.8 


32:03:11 


11.7 





16 


0.8 


13.8 


1.5 


1.7E+05 


100 


3:43:57.8 


32:04:06 


12.9 





16 


0.6 


13.3 


1.0 


1.6E+05 


105 


3:44:05.1 


32:00:28 


10.7 





20 


0.4 


8.2 


0.8 


9.1E+04 


109 


3:44:05.3 


32:02:05 


10.9 





14 


0.4 


16.0 


0.9 


1.7E+05 


110 


3:44:14.6 


31:57:59 


10.6 





12 


0.8 


11.1 


1.8 


6.9E+04 


111 


3:44:36.4 


31:58:40 


10.8 





15 


0.5 


10.5 


1.1 


1.1E+05 


115 


3:44:48.8 


32:00:29 


10.8 





10 


0.4 


9.4 


0.8 


1.6E+05 


117 


3:44:56.1 


32:00:32 


10.9 





12 


0.3 


6.8 


0.7 


7.4E+04 


118 


3:45:15.9 


32:04:49 


10.7 





08 


1.3 


11.6 


2.8 


4.9E+04 


119 


3:47:33.5 


32:50:55 


10.0 





08 


0.7 


9.9 


1.6 


5.0E+04 


121 


3:47:38.6 


32:52:19 


10.1 





17 


1.3 


18.5 


3.2 


1.8E+05 


122 


3:47:39.7 


32:53:57 


10.1 





10 


1.3 


18.5 


3.2 


1.8E+05 


122 


3:47:39.8 


32:53:34 


10.0 





11 


1.3 


18.8 


3.2 


1.8E+05 


122 



1 J2000 position of NH3 pointing, 2 J2000 position of NH3 pointing, 3 Kinetic temperature of core derived 
from NH3 (1,1) and (2,2) spectra, 4 Non-thermal 1-D velocity dispersion of NH3 spectra, °Total 1.1mm flux 
of the Bolocam core, Peak column density of the core derived from the peak 1.1mm flux and T kin , 7 Total 
mass of Bolocam core derived from the integrated 1.1mm flux and Tkin, 8 Averagc density of free particles of 
Bolocam core derived from its mass and size, 9 Per-Bolo number taken from Table 4 in Enoch ct al. (2008) 
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TABLE 2 

Environment of Starless Cores 



RA 1 


Dec 2 


^ km 


T d 4 


Too 5 


NIR A v 6 


FIR A v 7 


o-CO 8 


Distance Class O/I 9 


In a Cluster? 10 


J2000 


J2000 


K 


K 


K 


mag 


mag 


km s~ 1 


pc 




3:25:07.8 


on n a nn 

60:24:22 


9.2 


14.6 


11.6 


5.0 


2.6 


0.5 


1.51 


N 


• > or" nn t 

3:25:09.7 


on no F" o 

30:23:53 


9.2 


14.6 


11.2 


5.0 


2.6 


0.6 


1.50 


N 


3:25:10.1 


on A a A ~t 

30:44:41 


10.0 


15.2 


12.9 


5.1 


4.0 


0.6 


0.19 


N 


3:25:17.1 


on i o ft o 

30:18:53 


10.3 


12.9 


12.0 


5.0 


7.2 


0.7 


1.28 


N 


• i or" o o n 

3:25:26.9 


on n "i F" o 

30:21:53 


9.1 


12.9 


12.7 


3.7 


6.7 


0.5 


1.20 


N 


3:25:46.1 


30:44:11 


11.2 


15.5 


11.0 


5.5 


5.1 


0.7 


0.12 


N 


o o r" a o o 

3:25:48.8 


on a n n a 

30:42:24 


9.1 


14.3 


11.2 


4.8 


5.4 


0.7 


0.20 


N 


o o r" / \ / 1 

3:25:50.6 


on /in n n 

30:42:02 


9.0 


12.5 


10.6 


4.8 


10.4 


0.7 


0.24 


N 


3:z7:U2.1 


on i ft no 

30:15:08 


10.4 


13.5 


13.1 


4.8 


3.8 


0.6 


0.39 


N 


. > O T O O C\ 

3:27:28.9 


on "IF" n a 

30:15:04 


10.7 


14.9 


NaN 


9.5 


2.9 


NaN 


0.16 


N 


o it /inn 

3:27:40.0 


on i n i o 

30:12:13 


10.5 


17.4 


NaN 


8.8 


2.8 


NaN 


0.04 


N 


. > o t F" r" n 

3:27:55.9 


on n c "i o 

30:06:18 


9.1 


13.7 


NaN 


3.1 


3.1 


NaN 


0.17 


N 


o o o n r" r" 

3:28:05.5 


on n c "in 

30:06:19 


9.4 


13.7 


NaN 


2.9 


3.1 


NaN 


0.16 


N 


o no on a 

3:28:32.4 


O "1 n A A o 

31:04:43 


10.6 


12.8 


18.8 


5.4 


8.3 


1.1 


0.14 


Y 


• 1 o O O O A 

3:28:33.4 


on "in of" 

30:19:35 


10.7 


13.0 


NaN 


4.5 


6.0 


NaN 


0.90 


N 


o no on ft 

3:28:39.5 


31:18:35 


11.7 


12.3 


18.6 


7.8 


19.1 


0.8 


0.05 


Y 


o n O A O F" 

3:28:48.5 


O "1 1 / ■ C\ O 

31:16:03 


11.3 


16.0 


16.2 


8.5 


0.0 


0.9 


0.15 


Y 


o n o F" n n 

3:28:52.2 


o "i "i o n o 

31:18:08 


13.8 


15.0 


19.4 


9.2 


10.7 
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^ J2000 position of NH3 pointing, 2 J2000 position of NH3 pointing, 3 Kinetic temperature of core derived from NH3 (1,1) and 
(2,2) spectra, 4 Dust temperature of the Perseus molecular cloud derived from FIR emission, 5 Gas temperature of the Perseus 
molecular cloud derived from 12 CO (1-0) emission, 6 Column density of the Perseus molecular cloud derived from NIR extinction, 
Column density of the Perseus molecular cloud derived from FIR emission, 1-D velocity dispersion of 13 CO (1-0), 9 Distance 
to nearest Class 0/1 protostar, 10 Is the core within a cluster (NGC1333 or IC348)? 
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TABLE 3 
Tkin Correlations 



Property 


Location 


r value 1 


t value 2 
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1 coefficient of correlation (r), 2 t statistic: r/sqrt[(l — r 2 )/DOF], 3 % confi- 
dence of correlation, 4 degrees of freedom. 5 distancc to a Class 0/1 source, as 
identified in Enoch ct al. (2008b) and plotted in Figure [2] 




Fig. 1. — Histogram of T^ n . The filled (grey) histrogram is the temperature distribution of all 50 of the starless cores in Perseus for 
which we have well-fit NH3 spectra. The dashed (red) and solid (blue) histograms are the temperature distribution of the starless cores in 
the field and inside the clusters, respectively. 
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Fig. 2. — T]ji n plotted against properties of the enviroment around the starless cores: (top left) distance to nearest Class protostar, (top 
middle) distance to nearest Class I protostar, (top right) distance to nearest Class or Class I protostar (middle left) dust temperature, 
(center) gas temperature, (middle right) column density of the molecular cloud derived from NIR extinction, (bottom left) column density 
of the cloud derived from FIR e mission and ( bottom middle) 13 CO line width (1-D velocity dispersion). In this plot Class cores and 
Class I cores are identified as in lEnoch et atl Il2008bl) . Starless cores in NGC1333 and IC348 are plotted in black, and field starless cores 
are plotted in grey. The la errors in the fit of are plotted for each core. 
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Fig. 3. — Tfcin plotted against properties of the starless cores: (top left) non-thermal line width (1-D velocity dispersion) of NH3, (top 
middle) Bolocam Flux, (top right) peak Ay, (bottom left) core mass and (center) average volume density. Starless cores in NGC1333 and 
IC348 are plotted in black, and field starless cores are plotted in grey. The la errors in the fit of T^ n and non-thermal linewidth are 
plotted for each core. 



